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ABSTRACT 

We study the strong gravitational lensing properties of galaxy clusters obtained 
from N-body simulations with different kind of Dark Energy (DE). We consider both 
dynamical DE, due to a scalar field self-interacting through Ratra-Peebles (RP) or 
SUGRA potentials, and DE with constant negative w = p/ p = —1 (ACDM). We have 
12 high resolution lensing systems for each cosmological model with a mass greater 
than 5.0 x lO^^/i~^M0. Using a Ray Shooting technique we make a detailed analysis 
of the lensing properties of these clusters with particular attention to the number 
of arcs and their properties (magnification, length and width). We found that the 
number of giant arcs produced by galaxy clusters changes in a considerable way from 
ACDM models to Dynamical Dark Energy models with a RP or SUGRA potentials. 
These differences originate from the different epochs of cluster formation and from the 
non-linearity of the strong lensing effect. We suggest the Strong lensing is one of the 
best tool to discriminate among different kind of Dark Energy. 
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1 INTRODUCTION 

The mounting observational evidence for the existence of 
Dark Energy (DE), which probably accounts for ~ 70% 
of the critical density of the Universe (Perlmutter et al. 
(1999); Riess et al. (1998); Tegmark, Zaldarriaga, & Hamil- 
ton(2001); Netterfield et al. (2002); Pogosian, Bond, & Con- 
taldi (2003); Efstathiou et al.(2002); Percival et al. (2002); 
Spergel et al (2003)), rises a number of questions concerning 
galaxy formation. The nature of DE is suitably described by 
the parameter w — p/p, which characterizes its equation of 
state. The ACDM model (w = —1) was extensively stud- 
ied during the last decade. Recently much more attention 
was given to physically motivated models with variable w 
(Mainini, Maccio, & Bonometto (2QQ3a)), for which a num- 
ber of N-body simulations have been performed (Klypin et 
al 2003, KMMB03 hereafter, Dolag et al. 2003, Linder & 
Jenkins 2003, Maccio et al. 2004). One of the main results 
of KMMB03 was that dynamical DE halos are denser than 
those with the standard ACDM one. In this work we want 
to analyze the impact of this higher concentration on the 
strong lensing properties of the cluster size halos. 

Was first noted by Bartelmann et al. (1998) (for OCDM, 
SCDM and ACDM cosmology) that the predicted num- 
ber of giant arcs varies by orders of magnitude among dif- 
ferent cosmological models. The agreement between data 



and ACDM simulation was tested by many authors (see 
Meneghetti et al 2000, Dalai et al. 2003, Wambsganss et al. 
2004) but the situation is still unclear. A direct comparison 
of arcs statistic with observational data is out of the scope 
of this work, what we want to point out is the capability 
of Strong Lensing to discriminate between different kinds 
of Dark Energy (a similar paper but for a different choice 
of the dynamical DE parameters was recently submitted by 
Meneghetti et al. (2004)). 

Here, using a Ray Shooting technique, we make a lens- 
ing analysis of dark matter halos extracted from N-body 
simulations of cosmological models with varying w arising 
from physically motivated potentials which admit tracker 
solutions. In particular, we focus on the two most popu- 
lar variants of dynamical DE (Wetterich (1988); Ratra & 
Peebles (1988); Wetterich (1995)). The first model was pro- 
posed by Ratra & Peebles (1984, RP hereafter) and it yields 
a rather slow evolution of w. The second model (Brax & 
Martin (1999); Brax, Martin, & Riazuelo (2000); Brax & 
Martin (2000)) is based on potentials found in supergrav- 
ity (SUGRA) and it results in a much faster evolving w. 
Hence, RP and SUGRA potentials cover a large spectrum 
of evolving w. These potentials are written as 
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V{ct>) = exp(47rG0^) SUGRA. (2) 

Here A is an energy scale, currently set in the range 10^- 
10^'^ GeV, relevant for the physics of fundamental interac- 
tions. The potentials depend also on the exponent a. The 
parameters A and a define the DE density parameter fioB- 
However, we prefer to use A and Qde as independent pa- 
rameters. Figure 10 in Mainini et al. (2003b) gives examples 
of w evolution for RP and SUGRA models. 

The SUGRA model considered in this paper has A = 
10^ GeV this implies w — —0.85 at z = 0, but w drasti- 
cally changes with redshift: w ~ —0.4 at z = 5. The first RP 
model (RPi) has the same value for A of the SUGRA model. 
At redshift « = it has w = —0.5; then value of w gradu- 
ally changes with the redshift: at z = 5 it is close to —0.4. 
Although the w interval spanned by this RP model covers 
values significantly above -0.8 (not favored by observations), 
this case is still important both as a limiting reference case 
and to emphasize that models with constant w and models 
with variable w produce different results even if average val- 
ues of w arc not so different. For the second RP model (RP2) 
we have chosen A = 10~* GeV, in this case the value of the 
state parameter at redshift z = is the same of SUGRA: 
w{z = 0, A = 10"® GeV) = -0.85. This model is certainly 
better in agreement with CMB constrains but it loses most 
of its interest from a theoretical point of view: such a small 
value of A has not any clear connection with the physics 
of fundamental interactions and so it has exactly the same 
"fine tuning" problem of the ACDM model. 

We have normalized all the models in order to have 
today the same value of the rms density fluctuation on a 
scale of 8 /i~^Mpc , that has been chosen as crs = 0.8. 



2 N-BODY SIMULATIONS 

The Adaptive Refinement Tree code (ART; Kravtsov, 
Klypin & Khokhlov 1997) was used to run the simulations. 
The ART code starts with a uniform grid, which covers 
the whole computational box. This grid defines the lowest 
(zeroth) level of resolution of the simulation. The standard 
Particles- Mesh algorithms are used to compute density and 
gravitational potential on the zeroth-level mesh. The ART 
code reaches high force resolution by refining all high den- 
sity regions using an automated refinement algorithm. The 
refinements are recursive: the refined regions can also be 
refined, each subsequent refinement having half of the pre- 
vious level's cell size. This creates a hierarchy of refinement 
meshes of different resolution, size, and geometry covering 
regions of interest. Because each individual cubic cell can be 
refined, the shape of the refinement mesh can be arbitrary 
and match effectively the geometry of the region of interest. 

The criterion for refinement is the local density of par- 
ticles: if the number of particles in a mesh cell (as estimated 
by the Cloud-In-Cell method) exceeds the level rithrcsii! the 
cell is split ("refined") into 8 cells of the next refinement 
level. The refinement threshold depends on the refinement 
level. For the zero's level it is nthresh = 2. For the higher 
levels it is set to nthresh ~ 4. The code uses the expansion 
parameter a as the time variable. During the integration, 
spatial refinement is accompanied by temporal refinement. 
Namely, each level of refinement, I, is integrated with its own 



Model 


A 

(GeV) 


Box 

(/i-iMpc) 


Np 


(fe-iMe) 


(/i-ikpc) 


RPi 


10^ 


320 


512» 


2.03 X 10"' 


4.8 


RP2 


10-8 


320 


512^ 


2.03 X lO^" 


4.8 


SUGRA 


103 


320 


5123 


2.03 X 10^0 


4.8 


ACDM 





320 


512^ 


2.03 X 10^° 


4.8 



Table 1 . Parameters of simulations 



time step Ao; = Aoo/2', where Aao = 3 x 10"^ is the global 
time step of the zeroth refinement level. This variable time 
stepping is very important for accuracy of the results. As 
the force resolution increases, more steps are needed to inte- 
grate the trajectories accurately. Extensive tests of the code 
and comparisons with other numerical iV-body codes can be 
found in Kravtsov (1999) and Knebe et al. (2000). The code 
was modified to handle DE of different types (Mainini et al 
2003b & KMMB03). 

We performed a low resolution simulation for each 
model with the following parameters: box size: 320 ft-^Mpc, 
number of particles: 128^, force resolution: 9.2 ft-^kpc, all 
the simulations have the same initial random seed so at 2 = 
the clusters are more or less in the same positions. Than we 
selected the four massive clusters in the ACDM simulation 
and re-run them with a mass resolution 64 times higher. 
The same clusters are also re-run with the same resolution 
also in the RP and SUGRA models. At the end we have 12 
lensing systems (each cluster can be seen by three different 
orthogonal directions) for each cosmological model, with a 
meiss resolution of 2.03 x lO^^/i-^M© and a force resolution 
of 4.8 ft-^kpc. A complete list of simulation parameters is 
contained in table 1. 



3 LENSING SIMULATIONS 

In order to compute arc statistics for the models discussed 
above, we adopted a technique similar to the one originally 
proposed by Bartelmann & Weiss (1994). We center the clus- 
ter in a cube of 4 h~^Mpc side length and study three lenses, 
obtained by projecting the particle positions along the co- 
ordinate axes. This grants us a total of 12 lens planes per 
model that we treat as though being due to independent 
clusters, for our present purposes. 

We then divide the projected density field E by the 
critical surface mass density for lensing 

so obtaining the convergence k. Here c is the speed of light, 
G is the gravitational constant, while Dl, Ds, Dls are 
the angular-diameter distances between lens and observer, 
source and observer, lens and source, respectively. Once we 
set the lens and source redshift, the value of the angular 
diameter distance depends on the cosmological model. We 
detail this point in the next section. In Figure 1 we show 
the convergence map for one of the cluster, whose length 
scale size is 4 h~^Mpc. The deflection angle due to this 2D 
particle distribution, on a given point x on the lens plane 
reads: 
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yj\ 



(4) 



k(x) 



Here ijj is the position of the j-th particles and A'^ is the 
total number of particles. 

As direct summation requires a long time, we sped up 
the code by using a P"^M-like algorithm: the lens plane was 
divided into 128x128 cells and direct summation was ap- 
plied to particles belonging to the same cell of x and for its 
8 neighbor cells. Particles in other cells were then seen as a 
single particle in the cell baricenter, given the total mass of 
the particles inside the cell. 

The deflection angle diverges when the distance between 
a light ray and a particle is zero. To avoid this unwanted 
feature we introduce a softening parameter e in eq:(4); the 
value e is tuned on the resolution of the current simulation. 

We start to compute d{x) on a regular grid of 256x256 
test rays that covered the central quarter of the lens plane, 
then we propagate a bundle of 2048x2048 light rays and 
determine the deflection angle on each light ray by bicubic 
interpolation amongst the four nearest test rays (see sec- 
tion 3.2 for further discussion on the effects of the adopted 
resolution in the lens mapping). 

The relation between images and sources position is 
given by the lens equation: 



(5) 



Ahk{x) = = 5hk 



and the local properties of the lens mapping are then de- 
scribed by the Jacobian matrix of the lens equation, 

dyh „ doh 

7^ ~ Ohk 

axk oxk 

The shear components 71 and 72 are found from A^k 
through the standard relations: 



71 (^) 



72(2) 



'-[Au(x) + A2i{x)], 



(6) 



(7) 



(8) 



and the magnifications factor /i is given by the Jacobian 
determinant. 



[All {x)A22 {x) - A12 {x)A21 {x)] ^ . 



(9) 



det A 

Finally, the Jacobian determines the location of the 
critical curves Xc on the lens plane, which are defined by 
deXAixc) = 0. Because of the finite grid resolution, we can 
only approximately locate them by looking for pairs of adja- 
cent cells with opposite signs of det A. Then, the lens equa- 
tions 

yc = Xc - a{xc), (10) 
yield the corresponding caustics ijc, on the source plane. 



3.1 Sources deformation 

For statistical purposes one has to distribute and map a large 
number of sources. We are interested in arc properties and 
arcs form near caustics; so for numerical efficiency we have 
to distribute less sources in those part of the source plane 
that are far away from any caustics, and more sources close 
or inside the caustics. We follow the method introduced by 




Figure 1. Convergence map for one of the ACDM clusters. The 
side length is 4 /i~^Mpc. 



Miralda-Escude (1993) and later adapted to non-analytical 
models by Bartelmann & Weiss (1994). In the previous sec- 
tion we have obtained the deflection angles for the Xij (with 
i,j = l,...,2048) positions on the lens (or image) plane, using 
the lens equation (10) we can obtain the corresponding po- 
sitions on the source plane yij{xij), as usual we call this 
discrete transformation the mapping table. 

We model elliptical sources with axial ratios randomly 
drawn from the interval [0.5; 1] and area equal to that of a 
circle with radius Vs = 0.5". We first distribute sources on 
a coarse grid of 32x32, defined in the central quarter of the 
source plane covered by the light rays traced (due to con- 
vergence only a restricted part of the source plane can be 
reached by the light rays traced from the observer through 
the lens plane). From the mapping table we have obtained 
the magnification {/j.), if it changes by more than one (ab- 
solute value) between two sources, we place an additional 
source between both, in this way we increase the resolution 
by a factor 2 in each dimension. For the n-th iteration of 
source positions the criterium to add additional sources is 
that magnification changes by 2"~^ . We repeat this proce- 
dure four times to obtain the final list of source positions. 
To compensate for this artificial increase in the source num- 
ber density we assign a statistical weig ht of 2^-" to each 
image of a source placed during the n-th grid refinement. 
On average we have about 15000 sources for each lensing 
system. 



3.2 Arcs Analysis 

To find the images of an extended source, all images-plane 
positions x are checked if the corresponding entry in the map 
table j/lies within the source: i.e. for an elliptical source with 
axes a, b and centered in {yl; t/2) it is checked if: 
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Figure 2. Magnification vs. length/ widtii ratio for ACDM clus- 
ters. These two quantities are equal for an isothermal sphere lens. 



iyi^ + (ji^<i, (11) 

where (2/1,2/2) are the components of the vector y. 

Those points fulfiUing the previous equation are part of 
one of the source images and arc called image points. We 
then use a standard friends-of-friends algorithm to group 
together image points within connected regions, since they 
belong to the same image (the number of images of one 
source ranges from 1 to 5 for our clusters). 

We measure arc properties using a method based on 
Bartclmann & Weiss (1994). The arc area and magnification 
are found by summing the areas of tfic pixel falling into the 
image. Arc lengths are estimated by first finding the arc 
center, then finding the arc pixel farest from the centroid 
as well as the pixel farest from this pixel. The arc length 
is then given by the sum of the lengths of the two lines 
connecting these three points. The arc width is defined as 
the ratio between the arc area and the arc length. 

In Figure 2 we plot the relation between length/width 
ratio and magnification; we found a good agreement with 
previous results obtained by Dalai et al. (2003). The scatter 
in this relation is due to local fluctuations in the surface mass 
density, highly distorted images are also highly magnified, 
but the converse is not always true. 

Before proceeding with our analysis we have performed 
some tests on the resolution adopted in our ray shooting 
code; figure 3 shows the fraction of sources (number of 
sources divided by the total number) vs. their length/width 
for different values of the resolution of the lens mapping 
grid N^j. (results are for the ACDM model with zi = 0.3 
and Zs = 1.0). If the resolution of the lens mapping is not 
high enough the critical curves are too small compared to 
the source size and a spurious cut off in the number of arcs 
with L/W > 10 appears. This cut off is totally artificial and 
it vanishes for Nhr > 2048. As figure 3 shows the results are 
stable also for a higher value of Nhr (4096), then in order 



Figure 3. Fraction of sources mapped in an arc vs length/width 
ratio of the arc for different values of the resolution of the lensing 
mapping. All curves are for the ACDM model with 2; = 1.0 and 
Zs = 0.3. 

to have a good compromise between resolution and compu- 
tational time we have adopted Nhr = 2048 in the following. 



4 ARCS STATISTICS 

In this paper we aim to compare the lensing properties of 
a given cluster as it appears in different cosmological mod- 
els. There are throe main features that affect the number of 
giant arcs: the concentration of the halo, the total number 
of lensing systems at a given redshift and the value of the 
critical surface mass density (Ecr)- 

As predicted analytically by Bartelmann et al (2002) 
(for constant w models) and first noted in numerical Nbody 
simulations by KMMB03, and then confirmed by Dolag et 
al (2003) and Linder & Jenkins (2003), the concentration 
of dynamical DE halos is greater than the concentration of 
ACDM ones. Here we use the same definition of concentra- 
tion of KMMB03: the ratio of the radius at the overdensity 
of the ACDM model (103 times the critical density) to the 
characteristic ("core") radius of the NFW profile, (see how- 
ever KMMB03 for more details). A greater concentration 
increases the probability of forming giant arcs. In Figure 
4 we report the density profile of the same halo simulated 
in different cosmological models. The RPi halo is clearly 
denser and more concentrated than the ACDM halo with 
the SUGRA halo laying in between; the RP2 halo (which is 
not shown in this plot) has a concentration parameter close 
to the one of the ACDM model. 

The expected number of objects with a mass exceeding 
4 X 10^''/i~^Mo (in order to produce multiple images) at a 
given redshift (in this case z — 0.4) can be estimated using 
a Press & Schechter formalism (see Mainini et al 2003a). In 
dynamical DE, objects form earlier than in ACDM, so we 
have more lensing systems per Mpc^ at z = 0.3. This can 
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Figure 4. Density profile of the same halo simulated in different 
model. The red curve is for ACDM , the black one for SUGRA and 
the blue for RPi. The halo has a virial mass of 6 X lO^''/i~^M0 . 



Figure 6. Value of the critical surface mass density Ecr in dif- 
ferent cosmological models, for all the curves we choose a lens 
redshift = 0.4 



D, 




Figure 5. Number density of halos with a mass greater than 
4.0 X lO^^h-^MQ for ACDM (low solid curve), SUGRA (dotted 
curve) and the two RP models. The value of A is 10"^ GeV for 
both RPi and SUGRA and lO'^ GeV for RP2. 



be taken into account by multiplying the number of arcs by 
1.3, 1.21 and 1.12 in RPi, SUGRA and RP2 respectively. (In 
Figure 5 we report the evolution with redshift of the mass 
function for a mass threshold of 4.0 x IO^/i^^Mq). 

The evolution of the scale factor a with time also de- 
pends on the model. This implies that, at a given redshift 
z = 1/a — 1, the angular diameter distance Dad is model 
dependent; in fact its value is given by: 



Ho J a V ^m,0 



da 



(12) 



Here c is the speed of hght. Ho and Q,m,o are the present 
value of the Hubble constant and the matter density pa- 
rameter and n_D_B(a) gives the evolution of the DE density 
parameter with the expansion factor. To compute flnE{a) 
for RP and SUGRA models we have used the analytical ex- 
pression of Mainini et al (2003b). In Figure 6 we show the 
value of the critical surface mass density for the adopted 
cosmological models. The different values for Ecr mean that 
a ACDM halo yields more arcs than a dynamical DE halo, 
if they have the same surface mass density. The effect of 
the different values of the angular diameter distance tends 
therefore to reduce the number of arcs in dynamical DE 
models. 

The first result of our analysis is shown in Figure 7, 
where we plot the fraction of sources (number of sources 
divided by the total number) vs. their length/width ratio 
for a cluster (lens) redshift zi = 0.3, where lensing is most 
efficient for a source redshift of 1.0 (as shown later in figure 
11). As expected, the RPi model produces more distorted 
images, due to its more concentrated halos. The SUGRA and 
RP2 models are quite similar for L/W > 10 and they lay 
in between RPi and ACDM which, as expected, produces 
less highly distorted images than the dynamical DE models. 
We want to underline that part of the higher lensing signal 
due to the higher concentration of dark matter halos in such 
models is canceled by the increased Scr value. This effect is 
clearly illustrated in figure 8 where we have computed the 
arc statistics in a SUGRA model using the ACDM critical 
surface density. As expected, we have obtained a result for 
SUGRA that is closer to the RPi one. 

As pointed out by many authors (Wambsganss et al. 
2003, Dalai et al. 2003) the number of arcs that a clus- 
ter is able to produce is strongly related to the redshift of 
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Figure 7. Fraction of sources mapped in an arc vs length/width 
ratio of the arc. Upper curve is RPi (blue) than from top to 
bottom we have SUGRA (red), RP2 (cyan) and ACDM (black) 
results. 



Figure 9. Arcs counts for two different values of the sources 
rcdshift in the ACDM model. 
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Figure 8. Same plot of Figure 7 for the SUGRA model where 
the arc properties are computed using Ecr(ACDM ) (solid line) 
and Y.^r{SUGRA) (dashed hne). 



the sources (although the strength of this effect is not yet 
completely understood, different authors found different re- 
sults). In Figure 9 we plot the arc number counts for the 
ACDM model for Zs = 1 (dashed line) and Zs — 2 (solid 
line). As in previous work we found that the number of arcs 
increases if we increase the source redshifts. Figure 10 shows 
the same results of Figure 7 but for Zs = 2. 

As expected, the total number of arcs increases in all 
cosmological models. Again we have a sort of hierarchy of 




0.0001 



5 10 15 

Length/Width 



Figure 10. The same of Figure 7 for Zs = 2. With this value for 
the source redshift the number of arcs increases in all the models. 



results according to what is expected from the dynamical 
evolution of dark matter halos in the corresponding cosmo- 
logical model. 

Moreover, due to the lower difference in the value of 
Ec,- for this sources redshift, the four models are better sep- 
arated, expecially the SUGRA and RP2 ones. A difference 
between these two models is somewhat expected even if they 
have the same value of state parameter today {w = —0.84), 
this arises from the different evolution of w. in SUGRA it 
drastically changes with redshift (w — —0.4 a.t z — 5) when 
it is more constant in RP2 (w — —0.72 a.t z — 5). 

As final result in figure 11 we show the evolution with 
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Figure 11. Number of arcs vs the lens redshift for two different 
thresholds of the L/W ratio: 10 and 7.5. Solid line is for ACDM , 
dashed for RPi, dot-dashed for RP2 and dotted for SUGRA (2:^ = 
1.0). 

the lens redshift of the number of arcs for two different 
choices of the length/width ration: 10 and 7.5 (the redshift 
of the sources is Zs — 1.0). On average the the RPi model 
is always above the others, instead the different between 
SUGRA and ACDM is more or less constant at all redshift 
and the lensing signal decreases rapidly for z > 0.45. 

The RP2 model has a sort of double behavior: it is close 
to ACDM for z > 0.35 but it is more similar to SUGRA 
for z < 0.35, we think that this bimodality is due again 
to the evolution of the state parameter in this model expe- 
cially if compared to the SUGRA one: the ration wrp^/wsu 
decreases with redshift towards unity at 2; = 0, so it is less 
different from ACDM at high redshift in respect to SUGRA. 

The peaks in the lensing signal have slightly different 
positions in the different models. As argued by other authors 
(Torri et al 2004, Meneghetti et al. 2005) this could be due 
to time offset between merger events in different dark energy 
cosmologies. 



5 DISCUSSION AND CONCLUSIONS 

Models with dynamical DE are in an infant state. We do 
not know the nature of DE. Thus, the state parameter w{t) 
is still uncertain. In view of this functional indetermination, 
at first sight, it could seem that the situation is hopeless. 

In spite of that, we can outline some general trends 
that result from our analysis: in dynamical DE models, halos 
tend to collapse earlier than in a ACDM model with the 
same normalization at 2 = 0. As the result, halos are more 
concentrated and denser in their inner parts (KMMB03). 
Starting from this finding we have explored the consequences 
of this higher concentration, on strong lensing properties of 
dark matter halos, in SUGRA and RP cosmologies. 

We found that RPi halos (obtained assuming the clus- 



ter abundance of the power spectrum and a value for the 
energy scale A in the range suggested by the physics of fun- 
damental interactions) produce a higher number of arcs with 
a L/W > 10 if compared to the standard ACDM model. 
This model (RPi) is marginal consistent with observations 
and its purpose is mainly to illustrate the principal effect of 
a dynamical dark energy component on arcs statistic. 

The second model we analyzed based on RP potential 
(RP2) is more realistic from an observational point of view 
but less motivated by theoretical arguments. This model 
produce about 50% more arcs with L/W > 10 than the 
ACDM one for zi = 0.3 and Zs = 1 but it is marginally 
distinguishable from ACDM for lensing system at mod- 
erate high redshift [zi > 0.35, fig 11) or for high redshift 
sources/arcs [zi = 0.3 and Zs = 2). 

The SUGRA model is always in between the 
ACDM and the RPi models and it produces about 70-80% 
more arcs than ACDM . This difference is almost constant 
both changing the sources and the lens redshift and it tends 
to disappear for zi > 0.6 (for Zs — 1.0) where all the lensing 
systems considered in this paper (M; ~ 5 x 10^^ ) are un- 
able to produce highly distorted images (except in the test 
model RPi). We also noted that part of the stronger lensing 
signal due to the higher concentration of halos in dynamical 
DE models is partially canceled by geometrical effects that 
increase the critical surface density in such models (fig. 6 
and fig. 8). 

As final remark we would like to stress that arc statis- 
tic is a powerful tool to investigate the nature of the Dark 
Energy. The forthcoming observational surveys(i.e. CFHT 
Legacy Survey, SDSS and others) will improve the statis- 
tic of giant arcs on the sky (for example the RCS-2 Survey 
(Gladders et al. 2003) will cover an area of 830 deg^ and 
is expected to produce 50-100 new arcs). Such an observa- 
tional material will provide a discrimination between DE 
cosmologies possibly allowing to constrain the A scale of the 
SUGRA and RP potentials. 
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